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ABSTRACT 



We present [OIII]500.7nm and Ha+[NII] images and long-slit, high resolution 
echelle spectra in the same spectral regions of Sa2-237, a possible bipolar plane- 
^ ' tary nebula. The image shows a bipolar nebula of about 34"extent, with a narrow 

. waist, and showing strong point symmetry about the central object, indicating 

it's likely binary nature. The long slit spectra were taken over the long axis of the 
nebula, and show a distinct "eight" shaped pattern in the velocity-space plot, 
and a maximum projected outflow velocity of Vexp =106 km-s~^, both typical of 
expanding bipolar planetary nebulae. By model fitting the shape and spectrum 
of the nebula simultaneously, we derive the inclination of the long axis to be 
70°, and the maximum space velocity of expansion to be <308 km-s~^. Due to 
asymmetries in the velocities we adopt a new value for the system's heliocentric 
radial velocity of -30km. s~^. We use the IRAS & 21cm radio fluxes, the energy 
distribution, and the projected size of Sa2-237 to estimate it's distance to be 
2.1±0.37kpc. At this distance Sa2-237 has a luminosity of 340 Lq, a size of 
0.37pc, and - assuming constant expansion velocity - a nebular age of 624 years. 
The above radial velocity & distance place Sa2-237 in the disk of the Galaxy at 
z=255pc, albeit with somewhat peculiar kinematics. 

Subject headings: planetary nebulae: bipolar 



1. Introduction 

Most low mass stars spend some time 
at the end of their hves as Planetary Neb- 
ulae (PNe) before they die away as cool- 
ing white dwarfs. The processes by which 
these spherically symmetrical low mass 
stars become -often highly asymmetric- 
PNe are not well understood. Especially 
strong asymmetries exist in the sub-class 
of bipolar PNe. Candidates for the for- 
mation of bipolar PNe include binary cen- 
tral objects, fast rotating single stars with 
strong magnetic fields, and interacting 
wind models. Although, for the more ex- 
tremely coUimated nebulae, a binary sys- 
tem is becoming a generally accepted ne- 
cessity. A glance at Balick (1987), Schwarz 
et al. (1992) and Gorny et al. (1999) 
shows that a majority of PNe are signifi- 
cantly asymmetrical. A recent overview of 
asymmetrical PNe is given in Kastner et 
al. (2000). 

Bipolar PNe have been shown to have 
properties very different from those of the 
general population of PNe (Corradi & 
Schwarz, 1995). They form an interest- 
ing group combining shocks and photoion- 
ization, dust and gas, and very high out- 
flow velocities with interesting kinematics. 
They have giant dimensions (several are 
> Ipc), central stars with higher mass (as 
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in type I PNe), and surprisingly, evidence 
shows they have lower than typical lumi- 
nosities (Corradi & Schwarz, 1995). In 
some cases, the combination of high tem- 
peratures and excitation with low lumi- 
nosity and the presence of a variable Near 
Infra-Red (NIR) source indicates the pres- 
ence of a compact, hot binary component. 
In the case of M2-9, measurable expan- 
sion parallax (Schwarz et al. 1997) and 
internal movements (Doyle et al. 2000) 
have been observed. Several bipolar PNe 
possess point-symmetry and are, strictly 
speaking, not plane-symmetric as are the 
other bipolar PNe. The point-symmetry 
can be spectacularly clear (Schwarz 1993, 
Lopez et al. 1993) and has been suggested 
to be due to precession in the binary orbit 
of the central object (Schwarz et al. 1992, 
Livio 2000). 

The bipolar nebula Sa2-237 (PN GOll. 1+07.0, 
17Hn2' -15°45'13"J2000.) which is 

briefly presented in a poster paper by 
Masegosa et al. (1999), is discussed in 
some detail in this paper. Information on 
Sa2-237 is sparse, a SIMBAD search only 
resulted in IRAS fluxes, and a heliocentric 
systemic radial velocity measurement of - 
81km-s^^. The 18 references since 1983 are 
all to catalogs or papers discussing samples 
of PNe; apart from the Masegosa poster, 
no papers have been found about the ob- 
ject itself. The kinematic and morpho- 
logical analysis of Sa2-237 presented here 
will hopefully advance the understanding 
of the interesting characteristics of bipolar 
PNe. 

2. Observations & Data 

The spectra were secured with EMMI, 
the ESQ Multi Mode Instrument, mounted 
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at the NTT, during the night of 2000.12.04. 
Using the long sUt echelle mode with grat- 
ing 14 results in a resolution of 0.004nm 
per pixel on the CCD. A slit width of 
0'/8 on the sky gives a spectral resolving 
power of 70000 spatially projecting onto 
2.2 pixels. The CCD was a Tektronix 
TK2048EB chip with 24^ pixels each 
covering 0'.'28 on the sky along the slit of 
300"length. The CCD has a RON of 5e, 
a dark current of 1.7e pix~^-hr~^, and was 
used with a gain of 2.17e-ADU~^ in the 
unbinned mode. The order of the echelle 
grating was selected by using an Ha filter 
(ESQ No.596) passing also the [Nil] hues 
at 654.8 & 658.3nm, or an [OIII]500.7nm 
filter (ESO No. 589). The exposure times 
were 4800s (in two 2400s exposures) for 
the Ha spectrum, and 2400s for the 
[OIII] spectrum. For details on the NTT, 
EMMI and it's filters and gratings, see 
www.ls.eso.org/lasilla/Telescopes/NEWNTT/ 

The images were also taken with EMMI, 
but used in the RILD imaging mode. The 
same ESO filters used for the spectra were 
used to acquire the images, and the expo- 
sure times were 500s for the [OIII] image 
and 300s for Ha. 

The observations were reduced in the 
usual way using the IRAF package (see 
www.iraf.noao.edu) and the resulting im- 
ages and spectra are shown in Figures 1 & 
2 respectively. 

Our images show a bipolar nebula with 
a narrow waist extending over about 34", 
and point symmetry about the central ob- 
ject. The lobes are oriented more or less to 
the E and W, and are brighter in the north- 
ern part of the E lobe and the southern 
part of the W lobe. In the extreme NW of 
the W lobe there is a faint loop visible, in- 



dicated by an arrow in Figure 1, which may 

have an even fainter counterpart in the op- 
posite lobe. Deeper imaging may be able 
to resolve this, and also determine if there 
is fainter emission further away from the 
central object. Several other bipolars have 
had faint emission detected far from the 
brighter parts of the nebula: e.g. MyCnl8 
(Bryce et al. 1997) & M2-9 Kohoutek & 
Surdej 1980) therefore it is worth looking 
for this around Sa2-237. The object is ex- 
tended in 2MASS NIR images, with the 
E lobe being the brighter, and the central 
object is much brighter than the nebula. 

The position of the slit used to acquire 
the spectra in Figure 2 is over plotted on 
the [OIII] image. The same slit position 
was used to acquire the Ha+ [Nil] 658.4 
spectrum. The slit width of 0'.'8 is plotted 
to scale. 

Additional data for this paper has been 
extracted from the literature. This in- 
cludes the following fluxes and magnitudes 

for Sa2-237. 

IRAS: 12^ = 1.17 25/i = 6.01 60/x 
= 16.56 100/x = 8.29 (Jy); Persi et al. 
(1987): J=11.85 H=11.05 K=10.72 L=8.7; 
2MASS: J=13.06 H=12.57 K=11.83; Tylenda 
et al. (1991): B=16.08 V=15.50 & 
H/3=2.5-10-^^W-m-2; Condon & Kaplan 
(1998): 1.4GHz (21cm)=5.8±0.5 mJy 

The Persi et al. and IRAS results 
are nearly co-eval (1984-1985) while the 
2MASS magnitudes are from 1998; note 
the differences of up to 1.5 mag between 
the NIR measurements over the time inter- 
val of 14 yrs, indicating variability of the 
source. 

The reported systemic radial velocity of 
Sa2-237 from Beauheu et al. (1999) is - 
81km-s-^ 
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Fig. 1.— The [OIII]500.7 (left) and Ha+[NII]658.3 (right) images of Sa2-237. The position 
of the sht is indicated in the [OIII] image, and it's width of O'.'S is shown to scale. The arrow 
indicates the faint loop discussed in the text. Each image is 70"on a side, with N f and E <— . 



3. Analysis 

Radial velocities were extracted from 
the more intense regions of the spectra la- 
beled from a to h in Figure 2. 2-D gaussian 
fits were used to determine the velocities 
listed in Table 1. 

2Vmax is the projected expansion ve- 
locity between the fitted regions with the 
largest velocity difference. We are aware 
that there is very faint emission at higher 
velocities but the SN ratio is too low to use 
these data for velocity calculations. Also, 
the smearing out due to seeing dictates the 
use of the fitted peaks for velocity determi- 
nations. This Vjnax forms a lower limit to 
the true space expansion velocity since the 
projection has not been taken into account 
yet. 

Note the difference in 2Yrnax between 
the [Nil] lines (211km-s-^) and the [OIII] 



(205km-s-i) & Ha (206km-s-^) lines. An 
explanation lies in the behavior of the ma- 
terial and the excitation lines. [Nil] and 
other low excitation lines (e.g. [01], [SII]) 
form in the outer regions of these nebulae, 
while [OIII] is typically found closer to the 
central object. The velocity of the mate- 
rial tends to be directly proportional to the 
distance from the central object, hence the 
higher velocities in the [Nil] line as shown 
in Sa2-237. As far as we are aware, this 
is the first time that this effect has been 
observed in a bipolar nebula; it is common 
in elliptical PNe. 

A small intensity asymmetry between 
the blue and red shifted loops in the spec- 
tra indicates the likely possibility of a small 
misalignment of the slit. Such an offset 
would show up as an intensity asymmetry 
due to the point symmetric nature of the 
nebula. An estimated offset of (y/S from the 
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Fig. 2. — The 2-D spectra of Sa2-237. The spatial direction is vertical, the spectra run 
horizontally, and each sub-frame is 0.83nm by 58". a) is the [NII]654.8nm line, b) is Ha, c) 
[Olll]500.7nm, and d) [NII]658.3nm. The letters a-h in the Ha spectrum indicate the various 
parts of the spectrum where we have extracted velocities by gaussian fits. These features 
are identifiable in all four spectral lines. The extracted velocities are listed in Table 1. 
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Table 1: Measured line velocities in km-s ^ 
with respect to the rest A 



nm 



Line 


[OIII] 


[Nil] 


Ha 


[Nil] 


Part 


500.685 


654.806 


656.280 


658.339 


a 


54 


44 


43 


43 


b 


14 
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-60 


-54 


-48 


-55 


d 


-60 


-58 


-57 


-59 


e 


145 


153 


149 


152 


f 


134 


143 


142 


143 


g 


129 


124 


127 


126 


h 


56 


58 


56 


57 


2V 

^ * max 


205 


211 


206 


211 



central object is sufficient to have caused 
this effect. Our conclusions, however, are 
not affected by this small offset. 

We have derived the systemic velocity 
as the centre of the point symmetry of the 
eight shaped velocity-space plot. This fol- 
lows directly from our modeling and is an 
accurate estimate of the system radial ve- 
locity. By applying the appropriate correc- 
tion we find the heliocentric radial velocity 
of Sa2-237 to be -30km-s-^ The LSR ve- 
locity for Sa2-237 is -16km-s~^. Note that 
this is significantly different from the ve- 
locity reported by Beaulieu et al. 1998 & 
Durant et al 1998). For objects at 1=11° in 
the Galactic plane we expect a small posi- 
tive velocity of about VLS'R=-|-10km-s~^ at 
2kpc. Our VL5i(;=-16km-s~^, possibly in- 
dicates somewhat peculiar kinematics for 
Sa2-237 but we still place the object in the 
disk, albeit a bit far from the plane for a 
bipolar (type I) PN at 255pc. With the 
Beaulieu velocity it would have been im- 
possible to kinematically place Sa2-237 in 
the disk. 



The fluxes and magnitudes listed in Sec- 
tion 1 are used to create the energy distri- 
bution of Sa2-237 in Figure 3. The NIR 
values plotted are those from Persi et al. 
1987, since they are nearly co-eval with 
the IRAS data. 
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Fig. 3.- 
237. 



The energy distribution of Sa2- 



There is clearly a major contribution 
from the IRAS fluxes compared to the vis- 
ible and NIR. This is addressed further 
below. Integrating the fluxes over wave- 
length, and correcting by a factor of 1.5 
(Meyers et al. 1987) for the unobserved 
parts of the spectrum, yields L=77d^L0 for 
the luminosity of Sa2-237, where d is the 
distance to Sa2-237 in kiloparsecs. 

4. Model and discussion 

Using a simple model for the shape and 
kinematics of a generic bipolar nebula, in 
the manner of Solf & Ulrich (1985), we flt 
the morphology and spectrum of Sa2-237 
simultaneously. The model has the follow- 
ing properties: the velocity of expansion 
is proportional to the distance from the 
center of the material, there is symmetry 
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about the equatorial plane, and no time 
evolution. The nebula therefore expands 
in a self-similar way. The parameters are: 
the polar expansion velocity, Vp, the ratio 
of polar to equatorial velocity, R^, the in- 
chnation to the line of sight of the polar 
axis, i, and the age-distance"-^, Pad- The 
expansion velocity parametrized func- 
tion of latitude in the nebula is given by: 

V(0) = Ve+(Vp-Ve)sin^(|0|) 

where is the latitude, Vg & Vp are the 
equatorial and polar expansion velocities 
respectively, and 7 is the shape factor. For 
more details see the paper by Solf & Ulrich 
(1985) whose parametrization we have fol- 
lowed here. 

By fitting the shape of the image and 
the spectrum simultaneously, we derive 
the following parameters for Sa2-237: 
The inclination, i, to the line of sight is 
70°and the maximum expansion velocity 
is 260km-s-^. 

We estimate the distance to Sa2-237 
from the observational and literature data. 
Comparision with M2-9, a well studied PN 
(although some authors consider it to be a 
symbiotic nova; Sa2-237 may well be such 
a nove as well) which has a similar mor- 
phology as Sa2-237 and a well-determined 
distance of 640pc (Schwarz et al. 1997), 
and with Mira, gives the following results: 

Scaling the M2-9 distance of to that 
of Sa2-237 by the object's 1.4 GHz radio 
fluxes gives 1.66kpc, by luminosity gives 
2.68kpc and by size gives 2.11kpc. A simi- 
lar comparison with the IRAS flux of o Get, 
which has an accurate Hipparcos distance 
of 168pc, and assuming that the cool com- 



ponent in Sa2-237 is also a Mira, made 
plausible by the >1.5mag variability of 
Sa2-237, gives 1.96kpc. The average dis- 
tance to Sa2-237 based on these four esti- 
mates, is 2.1kpc, with a standard deviation 
of 0.37kpc, which can be taken as the for- 
mal error on this value. We therefore adopt 
a distance of 2.1kpc±0.37kpc to Sa2-237. 

The distance of 2.1kpc is used to com- 
pute the following parameters for Sa2-237. 
The size of the object is 0.37pc, it's lumi- 
nosity is 34OL0, and the age of the neb- 
ula, assuming constant expansion veloc- 
ity, is 624yrs. Using the model-derived 
inclination of i=70°, the space velocity of 
expansion of Sa2-237 is 308km-s~"^, using 
the maximum measured expansion veloc- 
ity 2V'""''. Note that this tends to slightly 
overestimate the space expansion velocity 
since the maximum projected velocity, for 
objects near the plane of the sky, is not the 
polar one, but more due to material at an 
inclination angle marginally smaller than 
the 70°derived from the model, and located 
at a latitude slightly below the poles of 
the nebula. This would bring the observed 
space velocity nearer that derived from the 
model (see also below). 

Gomparing the results derived from the 
observational data with the model param- 
eters, we get that the agedistance param- 
eter, Pad=0.30, in good agreement with 
Pad=0.28 from the model, given the un- 
certainties in all parameters. The com- 
puted expansion polar velocity at 1=70" is 
308km-s~^, cf. 260km-s~^ from the model 
fitting. Note that only 4° difference in i, 
that is 1=66°, gives 260km-s~^ instead of 
308km- s"^, making also this a reasonable 
agreement. We therefore adopt a space ex- 
pansion velocity of < 308km- s^^ for Sa2- 
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Fig. 4. — Model fits to the images and spectra of Sa2-237. Top: left the [OIII] image, right 
the Hq! image with model nebula over plotted; bottom: left the [OIII] spectrum, right the 
Ho; spectrum with model spectra over plotted. 
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37. 

All these parameters are typical of bipo- 
lar nebulae, and this strengthens our confi- 
dence that Sa2-237 is indeed a bipolar PN 
or symbiotic nova such as M2-9, He2-104, 
or BI Cru (Corradi & Schwarz (1993) and 
various others. 

The orientation near the plane of the 
sky of Sa2-237 implies that the central ob- 
ject is likely obscured by equatorial ma- 
terial, which would cause a large extinc- 
tion. There must be a pole to equator den- 
sity gradient to produce the bipolar shape, 
with the highest density at the equator. 
This is the reason why the IRAS contribu- 
tion is much larger than the optical-|-NIR 
flux, as the FIR radiation comes from the 
equatorial disk as reprocessed optical and 
NIR flux absorbed in the torus. The en- 
ergy distribution of M2-9 is similar and 
this object is also near the plane of the sky. 
Note that radio fluxes are not (or much 
less) affected by this and the distance ra- 
tios computed from the radio and IRAS 
fluxes conflrm this: the radio distance is 
smaller (1.66kpc) than the IRAS (~ lumi- 
nosity) distance (2.68kpc), indicating some 
extinction in the IRAS flux but none in the 
radio. 

One can, over the next years, attempt 
to measure the expansion parallax for Sa2- 
237, as we have done for M2-9, as it should 
be (y.'053 yr~^, at our distance of 2.1kpc. 
Taking one or two images per year would 
give a result within a few years, and would 
allow a direct check on our derived dis- 
tance. We have started such a program 
of observations. 

In the images, the nucleus of the neb- 
ula is relatively bright, but in the spec- 
tra no strong Hck emission appears. It is 



likely that it is a stellar spectrum, not an 
emission line core, as in symbiotic stars. It 
would be interesting to take a spectrum of 
this star and attempt to classify it. 

5. Summary & conclusion 

We summarize the properties of Sa2- 
237. 

It is a bipolar planetary nebula, simi- 
lar to M2-9, Ml-16 and others, projecting 
34"onto the sky, at a distance of 2.1kpc, 
with a size of 0.37pc, L=34OL0, and ex- 
panding at 2Vea;p=616km-s~^. It's point 
symmetric morphology hints at the pres- 
ence of precession and thus the central ob- 
ject may well be a binary. The object's 
inclination is 70°to the line of sight. 

Based on the symmetry in the observed 
expansion velocities, we adopt a helio- 
centric radial velocity for Sa2-237 of - 
30km-s~\ and place it in the Galactic disk 
at z=255pc. 

We suggest that the expansion parallax 
can be measured on a relatively short time 
scale as it is Ri(y/05 p. a. This would conflrm 
or adjust the distance to Sa2-237. 

Spectra of the central object may tell 
us what properties it has, as it seems to be 
stellar in nature and not an emission line 
core, as in many other such objects. 
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